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Abstract. In this work we study the Wigner functions, which are the quantum 
analogues of the classical phase space density, and show how a full rigorous semiclassical 
^ , scheme for all orders of h can be constructed for them without referring to the 

actual coordinate space wavefunctions from which the Wigner functions are typically 
calculated. We find such a picture by a careful analysis around the stationary points 
of the main quantization equation, and apply this approach to the harmonic oscillator 
^ ' solving it for all orders of h. 

in ■ 
o 
o 
^ , 

Q i 1. Introduction 

(N 

The Wigner functions (WFs) help us to picture the quantum states, that are typically 
represented as wavefunctions only in either configuration or momentum space, in the 
Ch ■ full phase space. They correspond to the classical phase space density. According to 
^ ■ the so-called Principle of Uniform Semiclassical Condensation (PUSC), they condense 
^ ■ on a classical invariant object (ergodic component) in the strict semiclassical limit 
■ h ^ 0, when they become predominantly positive on this effective support (Berry 
1977b, Robnik 1998). They are of extreme importance when trying to compare and 
relate the results of quantum mechanics to the classical ones. 

We typically obtain WFs by first finding the eigenstates in one of the usual 
representations from which we then calculate the WFs as e.g. in equation @. It is, 
however, an intriguing question whether it is possible to handle the WFs as independent 
objects in the phase space without referring to the corresponding eigenf unction. Such 
an approach was hinted at already by Heller (1976,1977). By a careful resummation of 
the Moyal bracket and a proper ansatz for the WF he managed to get an expression of 
its time evolution. Interestingly enough, this result does not reduce to a simple Liouville 
equation, the reason being the singular behaviour of the WFs in the strict semiclassical 
limit. On the other hand. Berry (1977a) calculated the semiclassical approximation to 
the WF by first using the semiclassical wavefunction, but since the end result can be 
expressed in a way that does not put either the coordinates or momenta into a privileged 



Semiclassical analysis of Wigner functions 



2 



position, these approximations to the WFs can be analyzed in the full phase space 
independently from the wavefunction approximations from which they were actually 
calculated. Ozorio de Almeida (1998) dealt with the Weyl representation in both 
classical and quantum mechanics, and managed to find a semiclassical periodic orbit 
formalism for the WFs that may be especially useful in the classically chaotic systems. 

Here we will try to find a quantum formalism that would expand the above ideas 
in a way that would enable us to deal with WFs completely independently from the 
eigenfunctions (in coordinate or momentum space), while at the same time we would 
like to expand their semiclassical picture to all orders of h. Osborn in Molzahn (1995) 
did a similar expansion for the Weyl symbols of operators, which are the generalizations 
of the WFs to operators other than the density operator. They, however, require that 
the symbols are regular in the semiclassical limit, which is not true for the WFs that 
have an essential singularity in this limit. Still, with their approach it is possible to find 
the phase space picture of the Heisenberg time evolution operator and act with that one 
on the irregular WF to get its time evolution. 



2. The Wigner-Weyl formalism 

We can represent the WFs within the broader Weyl formalism by which operators are 
assigned symbols that are functions of the phase space coordinates and momenta. The 
Weyl representation of an operator A is given by 

A{q,p) = / {q + x/2\ A\q — x/2) exp {—ipx/h) dx. (1) 



If A is self adjoint, the symbol A{q,p) is real. Also, by integrating over p and then q 
one can see that 

Tri = Jdq {q\ A\q) = ^ j dq dp A{q,p). (2) 

The WF by definition is just the Weyl symbol of the density operator p = lip) {ipl divided 
by 2TTh, 

W{q,p) = — - [ ip\q — x/2)ilj{q + x/2) exp {—ipx/h) dx. (3) 
It has a nice property that 



dq dp W{q,p) = 1, (4) 

which follows from (|^), meaning that the WF is properly normalized over the whole 
phase space. Normalization (^) follows also clearly from (|^) 

Operators can be represented as being elements of a linear space. We can find a 
basis and a scalar product within this space that will make the manipulations of Weyl 
symbols easier. We can see that the trace of the product of an operator with the adjoint 
of another operator, 

[A,b) = Tt{Ab^), (5) 
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indeed satisfies the conditions for it to be a scalar product of the two operators. This 
scalar product is real if the operators A and B are self adjoint. 
A proper basis for our work is the family of operators 

u}{q,p) = [ \q + x/2) {q — x/2\ exp {ipx / h) dx . (6) 



By taking into account {q\q') = S{q — q') one can show that these operators are self 
adjoint, meaning that 

{qi\i2}{q,p) 1^2) = y^*^ (^^Y^ ~ ^) 6^P(^P(^i ~ ^2)/^) = (fel ^{q,p) \qi))^ , (7) 
and therefore 

uj{q,py =u;{q,p). (8) 

Here 6{x) is the Dirac delta function. These operators are also orthonormal with respect 
to the chosen scalar product, 

{uj{qi,Pi),u;{q2,P2)) = 

{x \ qi+ Xi/2) {qi - Xi/2 \ + 2:2/2) (^2 - 3:2/2 | x) ■ 



2nh 

■ exp {i{piXi + p2X2)/h) dxidx2dx = 6{qi - q2)S{pi - ^2)- (9) 
In deriving this relationship we only have to use the property 

6{x — a)6{x — b)dx = 6{a — b) (10) 



of the Dirac delta function. 

With the help of the above expression the Weyl symbol of an operator A can be 
written as 

A{q,p) = V2^{A,Lj{q,p)). (11) 

Since the operators Cj{q,p) form a complete set of orthonormal operators, we can also 
write 

-4 = -^^^ J M(l^Ppi(l^P)dqdp, (12) 

which can be verified by insertion into ([ll|). This relationship is most helpful when one 
wants to find how the Weyl symbol of the product of two operators can be expressed 
by their respective Weyl symbols. Let 

C = AB. (13) 

The Weyl symbol of the operator C is therefore 



C(g3,P3) = V2nh{AB,uj{qs,P3)). (14) 
By substituting 

^ = y= J A{qi,pi)u}{qi,pi)dqidpi (15) 
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and 

B = — i= / B{q2,P2)u}{q2,P2)dq2dp2, (16) 
V /vrn. J 

we obtain 

C{q3,P3) = J A{qi,pi)B{q2,P2)TT{Lj{qi,pi)u{q2,P2)oj{q3,P3))dqidpidq2dp2. (17) 

After a rather straightforward derivation not unhke {^) we obtain 

C{q3,P3) = (^-^ J dqidpidq2dp2A{qi,pi)B{q2,P2) ■ 

■ exp {2i [pi(g3 - ga) + ^2(^1 - gs) + ^3(^2 - gi)] /^) • (18) 

The equation that determines the Weyl symbol of a product of two operators is therefore 
an integral equation which makes it nonlocal. This will be the main equation that will 
be dealt with in the following analysis of WFs. 

3. WKB expansion of Wigner functions 

We are now prepared to tackle the analysis of the WFs. We will be dealing with the 
stationary problem of quantum mechanics, which in the standard picture leads to the 
search for eigenenergies and eigenstates of the Hamiltonian operator. In this standard 
picture, the main equation which an eigenstate \ip) must satisisfy is 

Hm = E\i;). (19) 

When dealing with WFs the core object we refer to is the density operator which, for 
the case of a pure eigenstate, is written as 

p = 1^) m ■ (20) 

To ensure a proper solution, the quantization condition for the density operator in the 
nondegenerate case actually needs to satisfy a pair of equations (Curtright et al 1998) 

Hp = Ep (21) 

and 

pH = Ep. (22) 

If we transform these equations to the Weyl formalism using (|18]) we obtain the pair of 
equations 

(71^) / ^'?i^^'i'^^2'^^'2P(gi,Pi)^^(g2,P2) exp (±iAi23/^) = Ep{q3,p3), (23) 

where 

A123 = 2 [qi{p2 - Ps) + q2{P3 - Pi) + q3{pi - P2)] , (24) 

which corresponds to four times the area of a triangle spanned by the points (g„,p„), 
where n = 1 ... 3, in the phase space. 
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In a way similar to the usual WKB approach we can write the Weyl symbol of the 
density operator as 

p{qn,Pn) = exp{ian/h). (25) 

The above may seem like a contradiction with the requirement that the WFs need to be 
real. We will, however, see, that the above represents just a part of the total solution 
and when all the parts are taken together the final result can indeed be made real. As 
in all the cases that follow, the index n represents the evaluation of the proper function 
in the point {qn,Pn)- The equation (^31) then becomes 



TT 



dqidpidq2dp2exp{i(f)i23/h)H{q2,P2) = E, 



(26) 



where 



0123 = CTl - CTS ± A123. (27) 

The approach to give us the main order h solution to the above problem is the 
integration in the neighbourhood of the stationary points of the phase (pus- The 
equations that determine these points are 

d , 



dXn 



n23 



0; X e {q,p},ne {1,2}. 



{2i 



From this we determine the conditions for the stationary points {qf'\p^i^) and {q2'' ,P2 
where {q3'\pf^) is the point in which we wish to determine the WF, as being 

Brh /<^rr(0)\ 

0, 



dqi 

d(j) 

dpi 

(90 

dq2 
d4> 

dp2 



^ dp J 

(0) 

-Pi 





[pr - 












[gf - 










, (0)1 

+ P3 


=0, 





0, 



■ (0) 
ql 



(0) 



0, 



(29) 
(30) 
(31) 
(32) 



where a^^^ denotes the lowest order h contribution to a, as the basic stationary point 
analysis cannot reach any further. The brackets (• • .)i denote the function within them 
to be evaluated in the corresponding stationary point {qf'Kpf'^}, where it is obvious that 
the points 1 and 3 are the same. We can now shift our origin to the chosen stationary 
point, 

(33) 
(34) 



Ql 
Pi 



Ql 
Pi 



(0) 

Q3 , 
P3 5 



Q2 = q2 



(0) ^ 

Q^T^ 



(0) , 

P2=P2-P3 ± 



V 9p J, 

dq ). 



1 



(35) 
(36) 
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Rewriting equation (|26| ) into the new coordinates we obtain 



— ) / dqidpidq2dp2 exp (i [cti ± 2{qiP2 - q2Pi)] /K) H{q2,P2) = E, (37) 



where 



and 



(Ti = CTi - as 



qi 



. dp . 



Pi 



(38) 



H{q2,P2) = H{q2,P2). (39) 

For the quantities denoted by ~ , the index m naturally denotes evaluation in the 
corresponding point {qm,Pm)- 

The analysis has so far been focused on the leading order h contribution. We can 
use this leading order approximation to expand the analysis to all orders in h, with the 
leading order of this analysis being the same as above, and we may write 



(40) 



n=0 



We also perform a Taylor expansion to all orders in variables {qn,Pn) for all quantities 
in equation (P7[), 



~(n) 



j=0 



H2 = T.Hi'\ 



(41) 



(42) 



1=0 



where the indices (j) and (/) denote the order of the homogeneous polynomials of 
the expansion with respect to the corresponding coordinates {qn,Pn)- These shifted 
coordinates are, just like in the leading order analysis, obtained using the component 
af'^ that represents the leading, zero order contribution in the expansion of ai with 
respect to h, as given in equations ( p3D -(|36|). It will soon become apparent why such a 
choice of coordinates is proper. 
By also taking into account 



exp(xj 



oo k 



k=0 



the equation ([37| ) can be written in a fully expanded form 



dqidpidq2dp2exp{±2i{qiP2 - q2Pi)/^) 



ttTl 

oo 1 

y - 

Lk=0 



(43) 



(44) 



^ nj=0 



n~{n)(j) 



E 



An important relationship to be used in the following analysis is 



J dqdpexp{±2iqp/h)q^p'^ = nh (±i'- 



(45) 
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It can be obtained by noting that 

,".exp(±2WS)^(T4)"^^;^^5p!^ (46) 

and 

J dqexp{±.2ipq/h) = 7rh5{p) (47) 
holds, and therefore we obtain 

J dqdpexp{±2tqp/h)q^p'' = J dp Ut-\ ^Kp)- (48) 

By an m-fold per-partes integration of the above expression we obtain the desired result 



The equation (0) shows that, upon integration, the factor 5m,n in the expression 
(^) eliminates all the contributions in the multiple sum of the expression (|4^ ) for which 
the order / in the expansion of H2 with respect to the homogeneous polynomials does 
not match the polynomial order of the product formed by the various Taylor expansion 
terms of the phase ai. These product terms stem from the /c-th order in the expansion 
of the exponential function and the subsequent evaluation of the fc-th power of the series 
that represents the full expansion of ai. This leads to the condition 

^ = EJp, (49) 

where jp represent the Taylor orders of those terms o'l"''*^-'^ in the expansion of a^"'' that 
form the chosen fc-th order product term. 

The main goal of this semiclassical analysis is to sort the various contributions of 
the equation (^ij) in terms of the orders of h with which they contribute. We denote 
the order of h by which each term contributes to the total result by o. We again make 
use of the equation (H3). By carefully comparing it to the equation (HI) we may see 



that for each contribution to the multiple sum/product in equation ( PI) its appropriate 
order of Ti is given by 



^ p=i 



(50) 



By also taking into account (^), we obtain 



0= -k + Y^ [ip + Up] . (51) 
p=i 

or, equivalent ly, 

k 

0= -k + l + J2'^p- (52) 

p=i 

It is very important to note that for each ci-j"''^^-'''^ the inequality 

Jp + np>2 (53) 
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holds. We can show this by first noting that jp > 1 holds, which follows from the fact 
that due to the construction (stationary point) of a the zero order Taylor contribution 
is equal to 0. Since, however, we are basing our analysis by an expansion around the 
stationary point of the leading order. Hp = 0, of the expansion of a with respect to h, 
this means that for Up = 0, the linear, jp = 1, Taylor contribution is equal to as well. 
The above inequahty leads to 

k<o, (54) 

which can be shown to hold true by 

k 

0= -k + J2 [jp + "-p] >-k + 2k = k. (55) 
p=i 

This inequality tells us that the product terms that contribute with a given order a in 
the h expansion can never comprise a greater number of factors k than is the chosen h 
order a. We can also show that 

Jp < 2o, (56) 



holds, which can be seen by showing, from (pl|) 



k 

jp < "^[jp + np] = o + k < 2o. (57) 
p=i 

From this it follows that at a given order a of the h expansion the solutions can be sought 
locally as the order of the derivatives involved can never be higher than the order a. 
Another important inequality to consider is also 

Up < o, (58) 

which follows from 

k 

Up < Up = o + k — I < 0, (59) 
p=i 

by also noting k < I]p=iJp = ^- The relationship ( |55| ) tells us that only those terms 
can contribute to a given order o of the expansion of the equation (|4^) with 
respect to h for which the order n in the expansion of 5" does not exceed o. 

All the above expressions lead to an important result that for each order o of the 
expansion of the equation (|4^) over h there is always a finite number of terms involved. 
Even though the basic expansion could have been done with respect to any point in 
the phase space, using the stationary point (s) is the only choice which leads to the 
properties as given above. Using the above properties the system becomes at least in 
principle locally solvable since the equation for evaluating each order of the expansion 
of a with respect to h contains only finite order derivatives of the quantities involved. 

It is also important to observe that, at each order o in the h expansion, the term 
with the highest order n (= o) can only be linear {k = 1) and contains the first {j = 1) 
derivative of a[^K This means that the gradient of a[^'' is, for each order n in the 
expansion of 5"i over the powers of h, determined by all a["' ^ for which n' < n. 
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Using the above knowledge we may now try to rearrange the equation (0) and 
therefore ( p6[ ) with respect to the orders of h. As we determined above, for each order 
only a finite number of terms should contribute. A properly reordered form of the 
equation (|26| ) is therefore 

(ti^) /^^i^^i^^2C?P2exp(±2i(giP2 - g2Pi)/^) ■ (60) 

■ E E E ^ (^) E n i^nr^':-^^^'' = E, 

where we already dropped the terms that do not contribute upon integration due to the 
relationships (^9]) and (|5lD not being fulfilled for them. The sum over {np,jp}k,i,o is to 
be understood as a sum over all such combinations of indices Up in jp that, for given k, 
I in o, match the conditions (|49|) and (plj). 

With limiting the classical Hamiltonian to the form 

H{q,p) = Tip) + V{q), (61) 

which is by far the most common and for which the Weyl symbol becomes equal to the 
classical Hamiltonian, the ordering of terms with respect to the order of h becomes a 
little simpler since all the mixed derivative contributions evaluate to zero in this case. 
After integration only the terms where the products of derivatives of a*^"-' with respect to 
only q or only p, multiplied by the derivatives of H with respect to p or q, respectively, 
are preserved. In this case it becomes simpler to evaluate the integration in the equation 
(0) and the separations of contributions with respect to the order of h can be done in 
a semi-closed form. The equation (^) therefore becomes 

(^)2+(^) J dqidpidq2dp2exp{±2i{qiP2 - q2Pi)/fi)- (62) 



oo o 2o 1 . ■ 




E m^r^ 




2 {np,jp}k,i.o P=l ^ ^ 




where the powers of pi and qi that stem from the contributions of a have already been 
joined and taken in front of the product symbol. The brackets (. . .)j denote the function 
in them to be evaluated at one of the the appropriate points {qf'\pf'^} which are given 
by ( PU[]5^ ) that were obtained via the stationarity condition using the lowest order cr^°^ 



of the expansion of a with respect to the powers of h. The leading term in the h 
expansion, which is given simply by {H)2 is handled separately due to its somewhat 
different nature. 

Using the equation (HSl) we obtain 



1 1^^ 



o=l k=l 1=1 



E = {H)2 + E(^^)" E(-l)'^ E (±2 j (63) 
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Figure 1. Determining the lowest order h contribution to the phase of the WF. The 
closed curve represents the manifold (curve) of constant energy. The gradient Vcr^°-' is 
orthogonal to the vector a, with the the actual value of the phase tr'-^-' being given by 
the area A (see the main text for details). 



dp'- 



k 

1 E n 

2 {npJp}k,i.oP=^ 



+ 



k 

1 E n 

2 {npJp}k,i,oP=^ 




which is the main equation to be solved, and is nicely sorted with respect to the orders 

of n. 

The lowest, zero order h contribution is an expression that looks trivial at first, yet, 
however, it is quite involved, 

(64) 



as it gives the pair of equations 



1 fda^\ 

2 \ dp )^ 



1 (da^^ 
dq 



E. 



(65) 



We can use this pair of equations to determine the local gradient of the leading order 
phase contribution. It is determined by the chord that can be spanned between two 
points on the curve (manifold) of the constant energy, H{q,p) = E, for which (^1,^1) 
is the center. The size of the gradient is equal to the length of this chord, while at the 
same time the gradient is orthogonal to it (see figure |1|). 

Using this gradient we may also determine the actual value of the function a^^\ The 
easiest way to do this is to find all such chords that are parallel to the one corresponding 
to the chosen point (gi, Pi) and lie between this point and the energy surface. The centers 
of these chords form a path s in the phase space that starts on the energy surface and 
ends in the point (gi,pi). The change of phase along this path is given by 



a 



(0) 



(66) 
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where I denotes the length of the chord that corresponds to the given point s and dh is 
the component of ds that is perpendicular to the chord. The value of this integral is A, 
which gives exactly the area of the region between the chord around a chosen point and 
the curve of constant energy. This result is the same as obtained by Berry (1977a) where 
the phase of the WF was determined using the leading semidassical approximation for 
the WF. 

The above equations typically give a pair of solutions. By properly connecting these 
solutions at the caustics along with taking into account higher order corrections leads to 
the quantization conditions and subsequent determination of the semidassical energies 
(see Berry 1977a for details). 

Let us now consider the higher order h equations. All the terms that contribute in 
the linear order of h give the pair of equations 

1 f fdH\ /9(7«\ fdH' 
.dp J ^ \ dq ) ^ \dq 



(68) 



ih 



^2 





while the next order is already a more involved expression 




Equations that correspond to higher orders are quite similar, and they contain higher 
orders of derivatives of both H and cr*^"\ while at the same time higher order products 
of (7i"^ are involved. 

4. Harmonic oscillator 



As is almost customary in quantum mechanics, the test example for any new method 
is the harmonic oscillator. By properly scaling the coordinates the Hamiltonian can be 
written as 



Hiq,p) = 



p + q 



(69) 
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As the Hamiltonian is quadratic in both the momentum and the coordinate, only those 
terms of the equation ( p3| ) can feature in its analysis that contain at most the second 
order derivative of the Hamiltonian and, consequently, the phase a. At the same time 
the Hamiltonian is symmetric with respect to rotations around the phase space origin, 
and the same is true of the solutions 

cr(") = a("Hr) (70) 

which depend only on the distance 



V+P' (71) 

from the phase space origin. 

Apart from the lowest order in the expansion over the powers of h, the proper 
equations for all orders n > 1 in the h expansion for this system are given by 

-n". (,(.)_) (-^ + -_ j ( i;^ ) = 0, (72) 



which is obtained from the pair of equations ([63| ) by rewriting them in the radial 
coordinates, where the fact that the solution is symmetric with respect to rotations 

reduces both of these equations to the expression above. We also introduced 

1 

7(r)=(^-r2)'. (73) 



(75) 



For the lowest order solution we use the equation ( poD to obtain 

^ = 27(r). (74) 

The next order in the expansion of a over powers of h is obtained by the equation ([72|) , 
which gives 

da^ _ 1 /V(r) 1\ 
dr 2 \ 7(r) r J ' 

and, after integration, 

(T^^^ = i [ln(7(r)) + ln(r)] (76) 

In the figure |^ we show the semiclassical approximations to the WFs (dashed) for 
various eigenstates by using the two contributions above along with the exact solutions 
(full lines) 

WM{r) = \ exp (-^-Pl { ^) (77) 



where Ln represents the Laguerre polynomial of order A^. We used the relationship 
(Robnik 1998) 

dqdpW\q,p) = ^. (78) 
to normalize the WFs. 
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y "-x^ Re(r) 

K 



Figure 2. We show the complex plane of the variable r. With the thick line we depict 
the cut in the complex plane that needs to be taken out of the area of definition of the 
phase derivative, a' , in order for it to be defined uniquely. We also show the contour 
V that goes around the singularities in the points r — {0, ±k}, where k — ^/2E/lj, 
and from which the quantization condition is determined 



To find these approximate solutions as well as to perform furtlier analysis tlie 
solution in the whole complex plane needs to be carefully defined. Due to the 



singularities of da'-'^^ /dr and da''^^^ /dr in the points r = {0, ±y 2£'/u;} and the fractional 
power expressions in both of them we can only make these derivatives uniquely defined 
on the whole complex plane with the cut as shown in the figure |^. This cut, on the other 
hand, just coincides with the main domain on which we seek the solution. Therefore we 
obtain two contributions on this cut that are the limits of the expressions as obtained 
by the limit of approaching the cut from the upper or lower side, and therefore these 
expressions correspond to various sections of the path V. It is interesting to note that 
this cut is actually essential if we want the whole solution to be made real. Therefore 
we used the two branches when constructing the total solution, which are obtained 
by taking the positive and negative value of the square root in the definition of 7(r) 
which then makes the total result real. By using both contributions we may write the 
approximation to our Wigner function as 

W{r) = Acos + 6^ exp(-(T(^)), (79) 

where A is a real constant. Evaluating the integral of the equation (|75| ) therefore gives 

= / dqdp W'^{q,p) ^ ^ / dqdp exp(-2(T(^)), (80) 

where the value of the square of the trigonometric function was replaced by its average 
which can indeed be done in the semiclassical limit ^ — »• where this function is rapidly 
oscillating. In our case this gives 

-^ = ^" (81) 
7(r)r 



J dqdp exp(— 2cr^^-') = 2tt J r dr 



and therefore 

1 
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We still need to determine the phase shift 6, which is altered every time we encounter a 
singularity of da^^^ / dr when traversing the path V as shown in the figure 0. Although 
the expression (|76|) tells us that the weight of the logarithmic contribution (which are 
responsible for the phase shifts) when traversing the point is twice as strong as that 
at the other singular points, we also need to take into account that traversing the path 
V we only do a half of the full enclosure of this singular point. Upon encountering any 
singularity along the contour V we therefore need to shift the phase by — 7r/2. 

If at the same time we demand that the total phase upon the full traversal of the 
contour V must change by an integer multiple of 27r, namely 27rM, as the WF, which 
is exponentially dependent upon this phase, must be singlevalued, this leads to the 
quantization condition which will be given in full detail later. The difference is that we 
now only take into account the two lowest contributions of the expansion of the phase 
with respect to h, although this already gives the exact result for the eigenenergies in our 
example. For odd M it can be shown that we obtain semiclassical approximations for 
the WFs that are odd with respect to the reflection of r, which, however, contradicts the 
initial observation that the WFs must be invariant with respect to rotations around the 
phase space origin. For the even solutions (M = 2A^), on the other hand, we find that 
the phase shift in the expression (^) needs to be = — 7r/4 for r > if a^^'^^r = 0) = 
is chosen. This yields the explicit expression of equation ([79D, 



N 



W^ir) = \ ' sin yJ2N + 1 - y^- 

nh^'!Ty^2N + 1 - y'^ ^ 

-{2N + l)^TCcos(-j=^==]+^Y, r = Vhy (83) 



for < ?/ < \^2N + 1. This is exactly the result one would obtain by 
approximating the exact solution (|7^) using the large N approximation for the 
expression exp{—y^/2)L]\f{y'^) as found in (Szego 1959). 

This now completes our approximate treatment of the WFs for the harmonic 
oscillator, namely the two lowest orders, and now we turn to the exact treatment of 
the energy spectrum by considering all orders. By using a straightforward yet lengthy 
procedure of induction it is easy enough to show that the expression below, when inserted 
into ([72|) , gives the correct solution to the problem. 



dr (2E _ 



where aj,"-* are unknown rational coefficients, except for n = 0, 1, where they are fixed 
by (^P. 

We may now try to calculate the spectrum. We obtain it by taking a certain energy 
E in the above equations and then trying to find such a value E so that the WF 

p = exp (85) 



Semidassical analysis of Wigner functions 



15 




is singlcvalucd. This does not necessarily mean that the value of the phase a/h needs 
to be singlevalued, as it may change by an integer multiple of 2tt when traversing any 
closed path in the complex plane without changing the value of p after such a traversal. 

It is interesting to note that using the classical WKB method to determine the 
eigenfunctions we may not directly link the condition of singlevaluedness to the condition 
of the solutions being square integrable. Yet the singlevaluedness condition yields the 



Semiclassical analysis of Wigner functions 



16 



correct values for eigenenergies in systems that are exactly quantum solvable. Therefore 
the question can be posed whether the two conditions are equivalent or does this only 
apply to the solvable systems that usually possess some other special properties like 
solvability by the factorization method (Infeld and Hull 1951) or other (Cooper et 
al 1995, Robnik and Salasnich 1997a, 1997b, Robnik and Romanovski 2000a, 2000b, 
Romanovski and Robnik 2000). For most of them we may find the appropriate quantum 
canonical transformations (Lahiri, Ghosh and Kar 1998, Veble 2001). As we will see, 
using the singlevaluedness condition yields the proper solution in our case as well. 

Let us now choose the closed path V in the complex r plane as given in figure 
^ that encloses all the singularities of our problem. These are found in the points 
r = {0, ±i/ — }. The change of phase along this path. 



(86) 
4D as 

(87) 



Ay) = ^ / cr'dr, 
n Jp 

is given by the residuum of a' at infinity. We obtain it by rewriting equation 



2\p 



dr 



UJ- 



1 + 



?,n-l2E 



+ ... 



The leading term of such an asymptotic series is of the order r"^*^"*"^, with all the other 
terms comprising a higher negative power of r. As the residuum is given by the prefactor 
to the term containing r~^, the above expression can have a nonzero residuum only for 
= 0, 1. By taking into account equations (^) and (^), the evaluation of these residua 
therefore yields 

-2E 



Alp = 271 



(88) 



By specifying that the above change of phase needs to be an integer multiple of 27r we 
obtain the quantization condition for the energy 



E = huj 



M 1 
T ^ 2 



^9) 



where M is a nonnegative number. These solutions, however, also contain those that 
yield WFs that are odd with respect to reflection of the r coordinate. Since the proper 
solutions need to be invariant with respect to rotations around the phase space origin, 
only the even solutions are the proper ones. This leads to M = 2A^ and therefore 



E = huj 



(90) 



By constructing the full semiclassical WFs we therefore solved the problem of the 
harmonic oscillator to all orders of h without referring to the actual wavef unctions. 
A similar analysis for the infinite potential well (1-dim box potential) is in progress 
(Veble 2001,2002) 
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5. Summary and conclusion 

By devising a full semidassical analysis of WFs we managed to rewrite quantum 
mechanics, that is typically considered in either only the momentum or coordinate 
representation, into an independent full phase space formalism. We obtained the full 
semidassical equations to all orders of h for these functions. This enabled us to solve 
the problem of the harmonic oscillator as an example. 

It is easy enough to generalize the equations themselves to more than one dimension. 
The problems arise when trying to solve for the main order contribution, as the mere 
condition of the appropriate chords lying on the energy surface yields infinitely many 
solutions. We need to take other conditions such as the singlevaluedness of the WFs 
with respect to all traversals in the phase space into account, and this is far from trivial 
to implement. Most likely such a procedure, if it is found, can function well only in 
classically integrable systems, or possibly for the regular states in the mixed systems, as 
nonintegrability and the chaotic motion associated with it break the ordered structure 
of the classical phase space which is most hkely necessary for the generahzation of the 
above procedure to work. Finding the extension of the approach to more than one 
degree of freedom is therefore the main goal of the work to follow. 
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